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We theoretically investigate inelastic transport through anisotropic magnetic molecules weakly 
coupled to one ferromagnetic and one nonmagnetic lead. We find that the current is suppressed 
over wide voltage ranges due to spin blockade. In this system, spin blockade is associated with 
successive spin flips of the molecular spin and depends on the anisotropy energy barrier. This 
leads to the appearance of a window of bias voltages between the Coulomb blockade and spin 
blockade regimes where the current is large and to negative differential conductance. Remarkably, 
negative differential conductance is also present close to room temperature. Spin-blockade behavior 
is accompanied by super-Poissonian shot noise, like in nonmagnetic quantum dots. Finally, we show 
that the charge transmitted through the molecule between initial preparation in a certain spin state 
and infinite time strongly depends on the initial spin state in certain parameter ranges. Thus the 
molecule can act as a spin-charge converter, an effect potentially useful as a read-out mechanism 
for molecular spintronics. 

PACS numbers: 73.63.-b, 75.50.Xx, 85.65.+h, 73.23. Hk 



I. INTRODUCTION 

The remarkable miniaturization in semiconductor- 
based microelectronics in the past decades rapidly ap- 
proaches its natural limits. In this context, electronic 
transport through single molecules has received much 
attention, 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 Since individ- 
ual molecules are about two orders of magnitude smaller 
than chip features presently attainable, devices composed 
of single molecules would present an important break- 
through. Besides, transport through single molecules al- 
lows one to study nonequilibrium processes in systems 
dominated by quantum effects, but having a limited num- 
ber of relevant degrees of freedom. This is particularly 
interesting in the presence of vibrations and local mag- 
netic moments. Transport through single molecules is 
thus of interest for fundamental physics as well. 

It is very promising to combine molecular electronics 
with spintronics, i.e., the idea to employ spins to store 
and process information, in order to build more efficient 
memory devices on the nanometer scale^&iii The com- 
bination of these two strategies could be realized with 
magnetic molecules, i.e., molecules with a local magnetic 
moment^ - 

Experimental work in this field has focussed on the 
Zeeman splitting of the Coulomb-blockade (CB) peaks 
as well as the Kondo effect ^Siii More recent transport 
measurements for a Mni2 derivative have exhibited fine 
structure of the CB peaks originating from magnetic ex- 
citations and regions of complete current suppression due 
to spin selection rules pi 3 - These experiments have also 
led to a number of theoretical works piiiii 5 ^ 1 ^^ Fur- 
thermore, transport through quantum dots coupled to 
magnetic leads has been studied extensivelyi2ii22i2ii24 In 



particular, it has been shown for paramagnetic^ and 
ferromagnetic^ leads that an external magnetic field, 
which lifts the spin degeneracy of energy levels, leads to 
voltage ranges above the CB threshold where the current 
is strongly suppressed due to different tunneling rates for 
spin-up and spin-down electrons^ 5 - However, magnetic 
molecules coupled to magnetic leads have received little 
attention. 

An essential requirement for spintronics is the ability 
to effectively control and detect the spin. In a recent 
paper we have shown that magnetic anisotropy is cru- 
cial for slow spin relaxation in magnetic molecules and 
can lead to giant spin amplificationii 6 . If the molecule is 
prepared in a magnetic initial state, the current can be 
highly polarized for an exponentially long time interval. 
Thus the spin moment transmitted through the molecule 
depends on the initial orientation of the molecular spin 
and can be much larger than the initial molecular spin 
itselfii 6 - This allows one to effectively read out the spin 
of a molecule coupled to two nonmagnetic leads. The 
strong anisotropy is expected to arise from the interplay 
of the ligand field (e.g., in porphyrin complexes or Mni2 
derivatives^) and spin-orbit coupling as well as from the 
interaction with image dipoles. 

In the present paper we study the inelastic transport 
through an anisotropic magnetic molecule weakly cou- 
pled to one nonmagnetic and one ferromagnetic lead cm- 
ploying the rate-equation approachi£ii2*iiiiii& This con- 
figuration is motivated by the possibility to switch the 
molecule to a predetermined spin state, i.e., to write the 
spin, by applying a bias voltage alone, in zero magnetic 
fieldii 6 " We show that the proposed configuration leads to 
interesting physics beyond the effect of spin writing, in- 
cluding the occurrence of large negative differential con- 
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ductancc (NDC) at high temperatures. This effect is dis- 
tinct from the NDC found at low temperatures in the fine 
structure of differential-conductance peaks due to inelas- 
tic processes In the low-temperature transport the 
anisotropy leads to the appearance of a finite window of 
bias voltages for which the current is large, whereas it is 
strongly suppressed on either side due to CB and spin- 
blockade (SB), 13 i 15 i 24 i 25 i 26 i 27 i 28 i 29 respectively. By SB we 
mean the suppression of the current due to small single- 
electron tunneling rates, in our case due to density-of- 
states effects. In our system the large local spin leads to 
interesting modifications of SB, as discussed below. 

We also find that the charge transmitted through a 
molecule prepared in a particular spin state depends 
strongly on this initial state in certain parameter regimes. 
In fact the difference of the total transmitted charge be- 
tween preparation and time t — > oo is an exponential 
function of a certain energy difference over temperature. 
This effect is related to the giant spin amplification men- 
tioned above^£ but appears in the charge channel. 



II. MODEL AND METHODS 

Our results are obtained for a molecule weakly coupled 
to two metallic leads. Relaxation in the leads is assumed 
to be sufficiently fast so that their electron distributions 
can be described by equilibrium Fermi functions. We as- 
sume transport to be dominated by sequential tunneling 
through a single molecular level with onsitc energy e and 
local Coulomb repulsion U. The full Hamiltonian of the 
system reads H = H mo \ + ffi ea ds + H t , wherein 

flmol = (e-eV g )n+^n(n-l)- J s ■ S - K 2 (S Z ) 2 



and 



B(s z + S Z ) 



(1) 



describes the molecular degrees of freedom, -fficads = 
SQ=L,RSk C r e ak a lk C r a "ka represents the two leads a = 
L,R (left, right), and H t = E«=l,r 
t* a cracker) describes the tunneling. Here, the opera- 
tor c\ creates an electron with spin a on the molecule, 
eje-f + cjc| and s = J2aa' < 4- <T o-er'C (T '/2 are the cor- 
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responding number and spin operator, respectively. J 
denotes the exchange interaction between the electrons 
and the local spin S. We assume the anisotropy of S to 
be of easy-axis type, K 2 > 0. For simplicity we consider 
identical g factors for s and S. The external magnetic 
field is applied along the anisotropy axis and a factor 
gfj,B has been absorbed into B. Unless stated otherwise, 
the value of the gate voltage V g is chosen as zero. a^ kcr 
creates an electron in lead a with spin a, momentum k 
and energy e ak . 

The eigenstates of the unperturbed Hamiltonian H mo \ 
fall into sectors with n = 0,1,2 electrons^. Since 
[S z ot , .ffmoi] = 0, the eigenvalue to of S z ot is a good quan- 
tum number. For n = 0, 2 we obtain 



e(2, to) = 2(e - eV g ) + U - K 2 m 2 - Bm. (3) 



For n — 1 and —S + 1/2 < m < S — 1/2 there are two 
orthogonal states with energies 

e ± (l,m) = e-eV g -Bm+--K 2 (m 2 +^j±AE(m) (4) 

where AE(m) = [K 2 (K 2 - J)m 2 + (J/4) 2 (2S + l) 2 ] 1 / 2 . 
For the fully polarized states with n = 1 and m = ±(S + 
1/2) the upper (lower) sign applies if K 2 — J/2 is positive 
(negative). 

Second-order perturbation theory in the hopping ma- 
trix element t a produces a Fermi's Golden rule expres- 
sion for the transition rates between two many-particle 
states \n) and m) of the molccule^iS* 1 ^* 1 ^ 
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Here, denotes the density of states of electrons with 
spin a in lead a, v uc is the volume of the unit cell, / de- 
notes the Fermi function, and the e„ are the eigenenergies 
of H nm \ of many-particle eigenstates \n). 

The matrix elements C^ m are defined by C^ m = 
(n\c a \m). The occupation probabilities F n (i) of states 
|n) are obtained by solving the rate equations, 



dP n 
dt 



( pmR - 



P n Rr, 



0- 



(6) 



For the steady-state probabilities the left-hand side of 
this equation vanishes. In the following we assume 
e = 1 eV and U = 2 eV. Then the doubly occupied molec- 
ular states have exponentially small probabilities at the 
voltages and temperatures we consider. Also, the densi- 
ties of states satisfy D\/D\ = 1 and Df/Df = 0.01, i.e., 
lead L is nonmagnetic and lead R is a nearly half-metallic 
ferromagnet. For simplicity we assume = D^. The 
local spin is chosen as S = 2. 

The current through lead a reads^ 



/" = Te^K-7i m )P" 



(7) 



where the upper (lower) sign pertains to a — L (R), n n 
denotes the occupation number of the eigenstate \n) and 
Bf n _ n contains only those terms in the rates, Eq. JSJ), 
that involve lead a. A rate equation approach is also used 
to compute the current noise spectrum of the system, 



s (<LT(i) 81* {0)), 



(8) 



e(0, to) = —K 2 m 2 — Bm 



(2) 



where 8I a (i) = I a (t) — (I a ) denotes the current fluctua- 
tions in lead a and (I a ) the steady-state current^ 
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III. RESULTS 

The anisotropy of the local spin partially lifts the 
degeneracy of the molecular energy levels with respect 
to the magnetic quantum number to, where to denotes 
the eigenvalues of the z component of the total spin 
Stot^^ 1 ^ This leads to a splitting of peaks of the differ- 
ential conductance dl/dV at low temperatures, as shown 
in Fig.QJa) for finite and in Fig. ^b) for vanishing mag- 
netic induction Bm^ The complicated fine structure in 
the vicinity of the degeneracy point Vq arises from the 
anisotropy of the local spin and the exchange interac- 
tion of the electrons in the molecular orbital and the lo- 
cal spin. Each fine structure peak corresponds to an- 
other transition becoming available for single-electron 
tunneling >i2iiii& Note that dl/dV is asymmetric with 
respect to the bias voltage, i.e., several fine-structure 
peaks have a significantly different intensity when the 
bias changes sign. This results from the breaking of spin- 
rotation symmetry by the ferromagnetic lead together 
with the spin selection rules for the tunneling processes. 

The V-V g map shows three transport regimes: At low 
bias the current is thermally suppressed due to CB, ex- 
cept close to the degeneracy point. The electron number 
on the molecule is constant, i.e., n = for V g < Vq and 
n = 1 for V g > Vq, since all transitions between differ- 
ent molecular charge states are energetically forbidden. 
At large positive bias the conductivity of the molecule is 
high, since the electrons in the right lead have enough 
energy to overcome the energy barrier between the n = 
and n = 1 states. This is the conducting regime. At 
large negative bias the conductivity of the molecule is 
low. This current suppression is due to the SB mecha- 
nism explained in the following. 

Two main definitions of SB are used in the literature 
on transport through quantum dots. The original defini- 
tion refers to the phenomenon that transition probabil- 
ities for single-electron tunneling vanish between states 
corresponding to successive electron numbers if the to- 
tal spins differ by more than AS = l/2. 13 i 15 i 26 i 28 i 29 The 
other, more general definition refers to the situation that 
the tunneling rate for electrons of one spin direction is 
strongly suppressed relative to the other, e.g., due to fer- 
romagnetic leads^*21 or Zeeman splitting^ In this case 
the system can be stuck in a particular molecular many- 
body state because the rates for leaving this state are 
small. 

In our case the SB is related to the second 
mechanismjSiSiSI but the interaction between the elec- 
trons on the molecule and the local anisotropic spin 
leads to modifications: As soon as the bias is suffi- 
ciently high, all transitions between the n = and the 
n = 1 multiplets are energetically possible, as shown in 
Fig. ^c) . Selection rules for sequential tunneling require 
Am = 1/2. Both spin-up and spin-down electrons hop 
onto the molecule with equal rates, since the densities of 
states -D| and in the left (incoming) lead are equal 
for both spin directions. On the other hand, spin-up elec- 
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FIG. 1: (Color online) Two-dimensional plots of the differ- 
ential conductance dl/dV in the vicinity of one particular 
degeneracy point at low temperatures, (a) dl/dV for Zee- 
man splitting B = 0.05 eV. (b) dl/dV for vanishing magnetic 
field. The model parameters are J = = 5meV. When 
the magnetic field is switched on, the CB regime and the SB 
regime are separated by a finite window of bias voltage for 
which the conductance is high. Since it is low (governed by 
Df/Df) in the SB regime, NDC must occur where the con- 
ductance drops again, (c) Scheme of the molecular energy 
levels with all allowed transitions involving the states with 
n = and n = 1 electrons in an external magnetic field. In 
the SB regime all transitions are energetically possible if the 
bias voltage V is large enough. However, the rates for a spin- 
down electrons tunneling from the molecule into the right lead 
(curved arrows with slash) are strongly suppressed due to the 
small density of states. 



trons leave the molecule much faster due to the polariza- 
tion of the right (outgoing) lead. Electrons keep flowing 
through the molecule until a spin-down electron tunnels 
in. This spin-down electron can leave the molecule only 
with a very small tunneling rate due to the low den- 
sity of states D^. On the other hand, it can rapidly 
leave the molecule as a spin-wp electron if the local spin 
is simultaneously reduced by unity. However, the num- 
ber of possible spin flips is limited, depending on the 
initial spin state. Therefore, the molecule finally ends 
up in the singly charged state with minimal spin, i.e., 
S'tot = — S — 1/2. Further electron tunneling is blocked, 
since the left lead is energetically unreachable and the 
right lead has a low density of states for spin-down elec- 
trons, cf. Fig.^c). 
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FIG. 2: (Color online) (a) Current- voltage characteristics in 
the vicinity of the CB threshold assuming K2 = J = B — 
0.05 eV. The inset shows the current for T = 0.02 eV over a 
broader voltage range, (b) Energy level scheme for the n = 
and n = 1 spin multiplets as a function of the magnetic quan- 
tum number m. The bias voltage is just high enough to al- 
low the transition between the many-particle states with spin 
m = S and m — S + 1/2. The steady-state current is highly 
spin-polarized, since the tunneling of spin-down electrons is 
thermally suppressed. This regime corresponds to the plateau 
with enhanced current in (a). 



Another interesting feature shown in Fig. QJa) is the 
appearance of a finite window of bias voltages between 
the CB and SB regimes for which the conductivity of the 
molecule is high. This obviously leads to NDC when the 
SB regime with low conductivity is entered. Figure EJa) 
shows the current-voltage characteristics in the vicinity 
of the CB threshold at zero gate voltage for finite mag- 
netic induction. The external field tilts the molecular 
energy levels with respect to the magnetic quantum num- 
ber m due to the additional Zeeman energy, as sketched 
in Fig. Hfb). This removes the degeneracy of the spin 
multiplets. For small bias voltages and low temperatures 
the current is suppressed due to CB. But as soon as the 
transition from the ground state of the n = multiple! 
with m = S to the lowest-energy state of the n = 1 
multiplet with m = 5 + 1/2 becomes energetically al- 
lowed, cf . Fig. Hfb) , the current increases to the plateau 
shown in Fig. (2Ja) . These two levels are then equally 
occupied. Moreover, the current through the molecule 
is highly spin-polarized, since spin-up electrons tunnel 



rapidly through the molecule whereas tunneling of spin- 
down electrons is thermally suppressed. When the bias 
is further increased to allow tunneling also of spin-down 
electrons, the molecule in several steps goes over to the 
state with n = 1 and m = — S* — 1/2 and the current 
is strongly suppressed by the SB mechanism discussed 
above. The first required transition from n = and 
m = S to n = 1 and m = S — 1/2 has the highest en- 
ergy so that all following transitions become active at 
the same bias^ The width AV of the window of bias 
voltages with large current corresponds to twice the dif- 
ference of these two excitation energies and is given byl& 

eAV = 2[2SK 2 +B-AE(S+l/2)-AE(S-l/2)], (9) 

where AE(m) = [K 2 (K 2 - J)m 2 + (J/4) 2 (25 + l) 2 ] 1 / 2 . 
This holds as long as D^/D^ <C 1 and the Zeeman en- 
ergy B is not too small — for B — > we have to take into 
account that the states with n = and m = iS become 
degenerate ground states. The observation of the current 
plateau requires the temperature to be small compared 
to eAV, see Fig.^a). 

At high temperatures the current steps broaden so that 
the plateau of high current vanishes. Nevertheless, the 
current remains large at bias voltages close to the tran- 
sition point from CB to SB, as shown in the inset in 
Fig. Ela). For vanishing magnetic induction, B = 0, 
there is no enhanced current at low temperatures, as we 
will explain below. Nevertheless, at high temperatures 
a pronounced current maximum develops close to the 
transition from CB to SB, as shown in Fig. EI a). For 
anisotropy K 2 > 0, the peak grows and broadens with 
increasing temperature. Since deep in the SB regime the 
current remains small, a broad region of NDC develops. 
This is remarkable, since NDC, which has already been 
observed or predicted for many systems, is usually a low- 
temperature effect. The high-temperature NDC is a dis- 
tinct effect, since it appears even when there is no NDC 
at low temperatures. 

Before we explain the high-temperature NDC, we note 
that the temperature dependence of the I-V curves is 
qualitatively different for vanishing anisotropy barrier 
and large exchange interaction, K 2 <C kT <C J, as shown 
in Fig. Etb). Most striking is the fact that the maximum 
of the peaks stays constant for T — > 0. For this situation 
the energy levels are shown in Fig. ^c). It is obvious 
that only a single energy difference is relevant as long as 
the higher n = 1 quartet is not occupied. If eV/2 ex- 
actly equals this transition energy, the Fermi functions 
in the rates, Eq. JSJ, for these transitions all equal 1/2 
and are thus independent of temperature. Consequently, 
the steady-state probabilities and current are also inde- 
pendent of temperature at this bias. For lower bias all 
transitions from n = to n — 1 are thermally suppressed 
and for higher bias the relevant transitions become fully 
active (the Fermi functions approach unity) and we en- 
ter the SB regime discussed above. In both regimes the 
current decreases for T — > QM> In the SB regime this 
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FIG. 3: (Color online) Temperature dependence of the 
current-voltage characteristics at high temperatures for van- 
ishing magnetic field, exhibiting NDC. In (a) we assume 
K2 = J = 5meV, corresponding to the spin multiplets 
schematically shown in Figs.Qand|21 except that now B = 0. 
In (b) we assume a large exchange interaction J = 0.1 eV 
and vanishing magnetic anisotropy K2 = 0. (c) Energy level 
scheme for the n = and n = 1 multiplets as a function of 
the magnetic quantum number m for the case in (b). 



leads to NDC. It is remarkable that the SB — a quantum 
effect — leads to large NDC at room temperature. 

For the case of K2 > and, more importantly, small 
Fig- Eta), the current at the maximum increases with 
temperature, since transport through the higher n = 1 
quartet contributes more and more. For the same rea- 
son the current maximum shifts to larger bias voltages. 
On the other hand, the current maximum vanishes for 
T — > 0. In Fig.^a) we have assumed vanishing magnetic 
induction and there is no window of enhanced current at 
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FIG. 4: (Color online) Current noise spectrum. Shown are 
the normalized correlation functions Sll(ui) / Sll(u = 00) for 
three different bias voltages representing the three transport 
regimes: V = 1.8 V (CB), V = 2.2 V (SB), and V = 7.0 V 
(conducting regime). The zero-frequency limit gives the Fano 
factor F = S^w— »0)/2e|/| which tends to 3 in the SB regime, 



if we assume the one lead to be half-metallic, D s / D 
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low temperatures: At small bias, in the CB regime, the 
steady state has predominantly n = and m = —S due 
to the asymmetric tunneling rates into the right lead. As 
the bias is increased the system directly crosses over to 
the SB regime when the bias equals the excitation energy 
to the state with n = 1 and m = —S — 1/2. 

The shape of the NDC peaks in Fig. EJb) for vanish- 
ing anisotropy can be fitted by a function of the form 
l/cosh 2 x: Coming from the CB regime, the increase 
of the bias voltage energetically enables tunnelling pro- 
cesses with a rate proportional to the Fermi function 
f[(E—eV/2)/kT], where E denotes the excitation energy. 
However, the probability for the molecule being trapped 
in the SB state is also proportional to f[(E — eV/2)/kT]. 
Therefore, the current also contains a factor of the form 
(l-f[(E-eV/2)/kT])+C, where the first term accounts 
for the probability of the molecule not being trapped 
and the constant C describes the small tunneling prob- 
ability for an electron with the "wrong" spin. C is 
determined by the ratio of densities of states D^/D^ 
in the ferromagnetic lead. Altogether, we obtain / ~ 
l/cosh 2 [(£: - eV/2)/2kT] + f[(E - eV/2)/kT]C, which 
is dominated by the 1/ cosh 2 a; term for small D^/D^. 
This should be compared to Ref. 0, where the shape of 
the peaks in the differential conductance due the open- 
ing of additional transport channels is considered. In our 
case the peak occurs in the current. 

Note that this simple argument does not apply to the 
NDC peaks shown in Fig. |3Ja). For nonzero anisotropy 
Ki, the ground state is no longer degenerate and the 
various transitions between states of the n = and n = 1 
multiplets have different energies, which leads to a more 
complex line shape. 

Experimentally, the current suppression originating 
from SB could be distinguished from other mechanisms 
by investigating the shot noise of the system^ Fig- 
ure 0] shows the current noise spectrum for three dif- 
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ferent bias voltages representing the CB, SB, and con- 
ducting regimes. The latter is here reached by ap- 
plying a large negative bias, which overcomes the SB 
by allowing conduction through doubly occupied states. 
We here consider the normalized correlation function 
5'll(w)/»S'll(w— >oo) which gives information about tem- 
poral correlations of tunneling events for the left lead. 
According to Eq. (JSJ, is the Fourier transform of 

the current-current correlation function. The constant 
part, SliX"^ ~~ * comes from the autocorrelation con- 
tribution, which is a positive 5-function at vanishing time 
difference. In the conducting regime the noise spectrum 
exhibits a minimum at zero frequency, since the corre- 
sponding contribution to the correlation function is neg- 
ative. Whenever the molecule is doubly occupied, fur- 
ther electrons cannot tunnel in due to the Pauli principle. 
This is the usual "antibunching" effect for fermions. The 
typical frequency scale apparent in Fig.0]for the conduct- 
ing case is the inverse of the typical time the molecule 
remains doubly occupied, which is of the order of the 
typical tunneling time r = (2-K\t a \ 2 D L v uc /h)~ 1 m^ 

In the CB regime, Coulomb repulsion hinders elec- 
trons from entering the molecule and the Fano factor 
F = S(u! = 0)/2e|/|, which is defined as the ratio of the 
zero- frequency noise to the classical Schottky result^ 
is close to unity. In this case, single electrons tunnel 
through the system in rare, uncorrelated events, and 
quantum correlations are unimportant. 

In contrast, in the SB regime the noise is enhanced 
at lj = 0. This maximum results from a bunching of 
the charge carriers, which means that several electrons 
tend to tunnel through the molecule within a short time 
interval. On the other hand, the average waiting time be- 
tween such events is comparatively long. Results for the 
Fano factor reveal super-Poissonian shot noise, F > 1, 
in the SB regime. It reaches the value F x 3 for one 
half-metallic lead, Df/Df -> 0. This factor can be un- 
derstood by taking into consideration that the Fano fac- 
tor contains information about the charge of the current- 
carrying particles and the quantum correlations between 
them. In the SB regime, the molecule is in the singly oc- 
cupied state and has minimal spin for most of the time. 
Electrons from the left lead cannot hop onto the molecule 
until a spin-down electron is emitted into the right lead. 
The rate for this process is strongly suppressed. How- 
ever, if the spin-down electron does leave the molecule, 
the probability for further tunneling processes is high, 
since electrons of both spin directions may then tunnel 
into the molecule. A current is flowing until a spin-down 
electron occupies the molecule again, which leads to the 
SB state. The Fano factor can be obtained from the fol- 
lowing expression^ 



F=(N) 



(t 2 ) - (t) 2 (N 2 ) - (N)' 



(N) 



(10) 



Here N denotes the number of electrons tunneling 
through the molecule in one "bunch" and t the wait- 
ing time between such processes. Since the molecule is 



in the same SB state between the events, they are un- 
correlated, giving an exponential distribution of times t. 
Since the tunneling is dominated by the three states with 
n = 0, m = —S and n = 1, m = —S ± 1/2, one finds a 
probability of 1/2 N for N electrons in the bunch. This 
leads to (N) = 2 and {N 2 } = 6. Since (t 2 ) = 2{t) 2 for 
an exponential distribution, Eq. (|10J) yields F = 3 for 
the Fano factor. This result is of the same fundamental 
origin as the super-Poissonian shot noise found for quan- 
tum dots in the SB regime in Ref. |2^, which considers 
the case K 2 = J = (no local spin) and nonmagnetic 
leads. In this case the SB is induced by application of 
a magnetic fields We note that Eq. (|10fl also gives the 
observed value F — 1/2 for the conducting regime. 

Bulka et al& have investigated a similar system, con- 
sisting of a quantum dot connected to two electrodes, 
where the dot and one of the leads arc nonmagnetic 
while the second lead is ferromagnetic. In such an F- 
N-N junction, a large spin can be accumulated on the 
dot. This spin accumulation causes NDC and strongly 
enhances the current shot noise due to spin noise acti- 
vated in the NDC regime^ In contrast, super-Poissonian 
shot noise in transport through anisotropic magnetic 
molecules arises from an interplay of spin-flip processes 
involving a small local spin and spin selection rules. 

As noted above, giant spin amplification can occur in 
transport through anisotropic magnetic molecules cou- 
pled to nonmagnetic leads^& At bias voltages close to 
the CB threshold the total spin transmitted from one 
lead to the other can become exponentially large at low 
temperatures, if the molecule is prepared in a magnet- 
ically polarized initial state at time t = 0. If one lead 
is ferromagnetic the steady-state current is highly spin- 
polarized so that it is necessary to consider instead the 
excess transmitted spin. Analogously, one can define the 
excess transmitted charge: The molecule is prepared in a 
specific state |n) at time zero and then evolves according 
to the rate equations © . Since it approaches the steady 
state exponentially, the excess charge 



dt[i a {t)-{n\ 



(ii) 



is finite. The main observation is that AQ" can depend 
very strongly on the initial state \n) if the steady state 
shows SB. Practically it is much easier to measure the 
charge accumulation in the leads instead of the excess 
spinjSi since the accumulated charge is conserved (except 
for leakage currents) whereas the spin is not. Moreover, 
it should be easier to employ the excess charge for further 
data processing. 

Results for the excess transmitted charge as a function 
of bias voltage are shown in Fig. [Sfa), where an initial 
state with n = and spin m = S is assumed. The fine 
structure close to the CB threshold basically originates 
from the anisotropy of the local spin. Most striking is the 
exponential enhancement of AQ" at voltages above the 
CB threshold, where the steady-state current is already 
nonzero (but still small due to SB). When the bias is just 
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FIG. 5: (Color online) (a) Excess transmitted charge AQ L 
as a function of bias voltage in the vicinity of the first CB 
step for different polarizations of the right lead: D^/D^ = 
0.5, 0.3, 0.1, and 0.01. We assume J = 4meV, K 2 = lmeV, 
vanishing magnetic field and an initial state with n — and 
m — 2. The steady-state current for Df /D^ = 0.1 is also 
shown. The inset shows AQ L as a function of the anisotropy 
constant Ki, where the bias voltage V corresponds to the 
arithmetic mean of the first two excitation energies, (b) En- 
ergy level scheme of the n — and n = 1 multiplets. The 
molecule is prepared in the initial state with spin m = S. 



large enough to allow the two transitions from the n = 
state with extremal spin to the n = 1 state with extremal 
spin, m = ±S — > ±(S* + 1/2), the system initially pre- 
pared in a state with maximal spin, m = S + 1/2, for 
an exponentially long time only performs transitions be- 
tween the two extremal state connected by a solid line 
in the level scheme in Fig. |S[b). Thus a sizable — and 
completely spin-polarized — current is flowing until the 
molecule overcomes the anisotropy barrier by thermal 
activation. After that it rapidly relaxes towards the SB 
state, which shows a small steady-state current. On the 
other hand, if the molecule is prepared in the state with 
minimal spin, m = —S — 1/2, it already starts out in the 
SB state and the current is always small. Thus in the first 



case an exponentially large excess charge is accumulated 
in the leads in addition to the excess spin. In this regime 
the proposed setup functions as a spin- charge converter, 
i.e., spin information can be read through a measurement 
of the excess transmitted charge, a much easier task. 

On the other hand, if the bias is large enough to al- 
low further transitions involving other levels, the system 
rapidly relaxes back towards SB so that AQ" is of order 
e. Note that the excess transmitted charge increases for 
increasing anisotropy of the local spin, as can be seen 
from the inset in Fig. EI a). 



IV. SUMMARY AND CONCLUSIONS 

We have studied transport through a single anisotropic 
magnetic molecule coupled to one ferromagnetic and one 
nonmagnetic lead. Our main findings concern the NDC, 
super-Poissonian noise, and a strong spin dependence of 
transmitted charge. 

Two types of NDC occur. One appears at low temper- 
atures in an external magnetic field. Here the Coulomb 
blockade and spin blockade regimes are separated by a 
finite window of bias voltages for which the current is 
strongly enhanced. This is due to the interplay of bias 
voltage, magnetic field, and magnetic anisotropy, which 
allows only transitions between two specific molecular 
many-body states and prevents the spin flips necessary 
for spin blockade. The other NDC effect appears at high 
temperatures in the vicinity of the Coulomb blockade 
threshold. This is a distinct effect — it can occur even 
when there is no NDC at low temperatures. It can be 
viewed as a manifestation of quantum effects (the Pauli 
principle and spin selection rules) at room temperature. 

Spin blockade is accompanied by super-Poissonian shot 
noise, as found earlier for nonmagnetic quantum dots. 
Furthermore, the total charge transmitted through a 
molecule under a voltage bias depends strongly on its 
initial spin state in a certain parameter range. The differ- 
ence in transmitted charge can diverge exponentially for 
low temperatures. This spin-charge conversion presents 
a promising method to read out the spin information in 
molecular-memory applications. 
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